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DEFINITION OF SYMBOLS 


Symbol 
A-l, A 

■^•oo » Am 

A(?,n) 

Bl 

B(5,n) 

c ml 

C(?,n) 

D(5,n) 

Em 

fo y El > f 2 
F 

G 00> 

g 

H 

Hi 

h 

Ejnl 

Jl 

k 

1 ,m 

L 2 

M 

M c 


Definition 

amplitude of incident and reflected 
sound pressure waves respectively 

sound pressure coefficients defined by 
Eqs. (39) and (47) respectively 

function defined by Eq. (68) 

sound pressure coefficient defined by 
Eq. (47) 

function defined by Eq. (69) 

sound pressure coefficient defined by 
Eq. (47) 

function defined by Eq. (73a) 

function defined by Eq. (73b) 

sound pressure coefficient defined by 
Eq. (51) 7 

non-dimensionalized perturbed sound 
pressures defined by Eq. (28) 

sound pressure defined by Eq. (20) 

sound pressure coefficients defined by 
Eqs. (41) 6(42) respectively 

transverse acoustic velocity defined 
by Eq. (22) 

duct height defined in Fig. (4) 

sound pressure coefficient defined by 
Eq. (44) 

vertical acoustic velocity defined by 
Eq. (21) 

sound pressure coefficient defined by 
Eq. (46) 7 

sound pressure coefficient defined by 
Eq. (45) 7 

sound wave number (=w/c 0 ) 

running integers 

parameter defined by Eq. (9) 

mean flow Mach number 

mean flow centerline Mach line 

sound field group Mach number defined 
by Eq. (65) 



Symbol 

N 


P 

R 


u,v,w 



x,y,z 

3 

5 n 


A 

K 

P 


Definition 

exponent on mean flow velocity profile, 
defined by Eq. (31) 

acoustic pressure 

amplitude of upstream sound pressure, 
defined by Eq. (7b) 

acoustic velocites in x,y,z directions 
respectively 

mean flow speed 

mean flow center line speed 

sound phase speed 

sound group speed defined by Eq. (63) 

duct width, defined in Fig. (4) 

duct coordinates, defined in Fig. (4) 

duct aspect ratio (sH/W) 

duct nondimensionalized transverse 
coordinates defined by Eq. (26) 

sound wavelength 

system eigenvalue defined by Eq. (20) 
density 


u sound circular frequency 

$ phase angle shift between incident and 

reflected sound waves in duct, defined 
by Eq. (7b) 


Subscripts 

o 

c 

max 


denotes mean quantities 
denotes duct centerline values 
maximum value of quantity 
minimum value of quantity 


mm 



SUMMARY 


The influence of a mean vortical flow on the connection 
between the standing wave pattern in a rectangular three-dimen- 
sional waveguide and the corresponding duct axial impedance has 
been determined analytically. The solution was derived using a 
perturbation scheme valid for low mean flow Mach numbers and 
plane-wave sound frequencies. The results show that deviations 
of the standing wave pattern due to refraction by the mean flow 
gradients are small. 
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INTRODUCTION 


Three different methods to measure the impedance of sound 
absorbing liners in the presence of grazing flow have been report- 
ed in the literature* 8 . The three methods are (1) the two-micro- 
phone method of measurement (References 1-6), (2) the impedance 
tube method (Reference 7) and (3) the wave guide method (Refer- 
ence 8). The two -microphone method is a particularly attractive 
method of measuring the impedance of cavity-backed, thin face- 
sheet liners consisting of either . porous materials or perforates. 
Its use however, is restricted to single cavity liners - thus its 
application is restricted to the measurement of the impedance of 
only locally reacting liners. 

The impedance tube method has been used successfully to 
measure the normal impedance of locally reacting materials in a 
nongrazing flow environment. In an attempt to apply this method 
in a grazing flow environment, Feder and Dean 7 experienced con- 
siderable difficulty in properly accounting for the radiation 
impedance losses of the sample liners mounted in the side of the 
flow duct facility. 

Armstrong 8 outlined a waveguide method to measure the 
impedance of a liner installed along one side wall of a constant 
area rectangular duct facility with grazing flow. By assuming 
that the rectangular duct consisted of a source at one end and a 
non -reflective termination at the other, Armstrong- was able to 
model the duct mathematically as semi- inf inite with one side 

lined, the others acoustically rigid. With this assumption, 
he proceeded to solve for the sound field within the duct assuming 
a uniform mean duct flow. In this way, Armstrong was able to 
show a connection between the liner impedance and the sound field 
attenuation and phase change per unit duct length (recall that 
one side wall of an entire infinite duct was lined - thus the 
sound pressure decayed continuously). Straight forward measure - 
of the sound field attenuation and phase change per unit 
length with an axially mounted microphone yielded, after relatively 
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minor computations, prediction of the liner impedance. Although 
this approach appears to be straight forward, difficulties 
associated with flow noise in the downstream diffuser (the flow 
appeared to separate in the diffuser) prevented Armstrong from 
assessing the accuracy of the waveguide approach. 

Despite the difficulties experienced by Armstrong, the 
wave guide approach is attractive because it will permit, at 
least conceptually, the measurement of the effects of flow on 
the impedance of extended liners. As noted above, the other 
two methods described above are restricted to locally reacting 
liners (by "locally reacting" from a practical viewpoint, it 
is meant that the size of the liner is very small compared to 
the wavelength) . A critical aspect in the development of a 
waveguide method of measuring impedance in a flow duct environ- 
ment is the connection between the standing wave pattern along 
the wave guide and the corresponding duct axial impedance. Of 
special interest here is the details of this connection when 
the mean flow is vortical. An analytical solution has been 
derived by restricting the Mach number of the mean flow to be 
small and the sound frequencies of interest to be below the 
first rigid-wall, zero Mach number cut-on mode - the plane wave 
mode. Under these restrictions, deviations of the standing 
wave pattern from its nongrazing flow plane wave distribution 
not be large. This permits an analytic representation 
of the flow related standing wave pattern in terms of regular 
perturbations from its zero grazing flow distribution. 

The use of the perturbation theory has been successfully 
used by Hersh and Catton 9 to study the effects of mean flow 
velocity gradients on the refraction of sound in infinite two- 
dimensional rigid-walled ducts. Their results showed good 
agreement with exact numerical solutions. 

Within the past decade, there has been an enormous interest 
in understanding the manner in which sound attenuation in lined 
ducts is affected by flow. A rather extensive review of this 
subject was recently presented by Nayfey, Kaiser and Telionus 10 . 



Briefly, previous studies in this field can be divided into two 
parts, one describing the effects of flow on the propagation 
of sound in rigid-walled flow ducts and the other to the propa- 
gation (and attenuation) of sound in acoustically lined flow 
ducts . 

Pridmore-Brown 11 was one of the first to study the effects 
of vortical flow on the propagation of sound in two-dimensional 
ducts. Since the publication of Pridmore-Brown ' s work, there 
has been considerable interest in the manner in which sound propaga- 
tes in a vortical flow. Recent publications of interest include 
the work of Mungur and Gladwell 12 , Mechel, Mertens and Schilz 13 , 

Tack and Lambert 14 and Mungur and Plumblee 15 . These studies 
which are primarily limited to the case of downstream propagation, 
have results consistent with those of Pridmore-Brown. 

It is intuitively clear that if the effect of mean flow 
shear is to refract thereby amplify the sound pressure near the 
duct wall for downstream sound propagation and into the duct cen- 
ter for upstream sound propagation, then the standing wave pattern 
and the corresponding duct axial impedance must take this into 
account . 

The report is organized as follows. The solution of the 
standing wave pattern for the case of uniform mean flow is 
presented in Section 2. The basic equations governing the 
propagation of sound in a rigid-wall three dimensional rectangular 
duct containing a mean vortical flow are derived in Section 3. 

The analytical solution is based on application of the perturbation 
method described in Section 4. A detailed solution is derived in 
Section 5 corresponding to the special case of plane wave propaga- 
tion in a square duct containing a fully-developed mean flow. 

The velocity profile of the mean flow was assumed to be linear. 

This permitted an analytic solution to be derived. The main 
findings of the study are summarized in Section 6. 
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UNIFORM FLOW SOLUTION 


The connection between the standing wave pattern in a wave- 
guide and the corresponding duct axial impedance is derived below 
for the special case of uniform duct flow. This case is important 
because it provides a check case for the more complicated vortical 
flow solution described in Section 4 and serves as an excellent 
introduction to the complexities associated with a vortical mean 
flow. The solution is presented in two parts, one in which the 
basic equations and the other the duct axial impedance are derived. 


2 . 1 Basic Equations 


Assume that the duct shown in Fig. 1 contains a uniform 
mean flow. Assume further that sound, introduced into the duct 
at, say, x=-L, generates a standing wave pattern between x=-L 
and x =0 . The strength of the reflected wave is a function of the 
duct effective axial impedance Z. The objective of the analysis 
is to derive a relationship between the standing wave pattern and 
the duct axial impedance Z with the duct mean flow speed, V 0 , a 
parameter. 

For uniform flow in a rectangular duct, the linearized con- 
tinuity and momentum equations are 


bp , \/ bp „ fb\k bv b w 

~T~ + *o — r~ + Po ( 4 1 

<H <}x ' \bx by 


du 


+ -L 

o 

n 

dt 

bx 

e° 

dx 

<dv 

+ V 0 ^l 

l 

+ — 

n 

o 

at 

* <bx 

Po 

a y 

b w 

+ 

+ — 

0 
II 

4 

1 


dx 


, b % 


( 1 ) 

(2a) 


(2b) 


(2c) 
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Here the subscript (o) denotes mean flow quantities. Thus p 0 and 
V 0 represent the mean density and uniform duct speed respectively. 
The quantities p,p,u,v,w represent acoustic excitations of pressure, 
density, axial, transverse in y direction, and transverse in z 
direction velocities respectively. Assuming p and p are adiabati- 
cally related, Eqs. (1) and (2) may be combined to form the well 
known convected wave equation 





ii 


& 

dz 2 


_ O 

C Q dx. dt 


(3) 


where M 0 =V 0 /co. 

For frequencies below the first cut-on mode, only plane 
wave modes are excited in a hard-walled duct. The acoustic pressure 
and velocity fields excited by the source are obtained by solving 
Eq. (3) to yield 



i (cot ~ k+x 

0 i (cot + k x) 


p ( x ' t ) 

= A +• e ' 

y +A-e v y 

(4) 

u (x, t) 

A + e A. e '(“ t + k - x ) 

(5) 

o 

o 

^0 

o 

n 

0 



where A+ and A_ represent the (complex) amplitudes of the incident 
and reflected waves respectively. The wave numbers 

= — — > ^ - = — — — ) K- = — — (6a,b,c) 

! + M 0 1-M 0 c 0 

represent respectively the elongation and contraction of the sound 
wavelength A = 2tt c 0 Ao due to the mean flow. Since the objective 
of this analysis is to derive an analytic expression connecting 
the impedance of a specimen located effectively at x =0 to the duct 
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mean flow speed, no loss of generality results by setting 


A + = I and A_ = Re 


(7a, b) 


The time-averaged pressure standing wave distribution along 
the duct follows immediately from Eq. (4)upon substitution of Eqs. 

(6) § (7), 

p 2 (x) = I + R. z + 2 R cos 

Equation (8) describes the output of a microphone at location 
x - see Figure 2. The impedance of the specimen is assumed to be 
located at x=0. It contains two unknowns, R the amplitude of the 
reflected wave and <j> the phase shift between the incident and 
reflected wave. Referring to the standing wave pattern shown in 
Fig. 2, L 2 is defined as 


2K; 


I- M : 


+ $ 


( 8 ) 


L z 


2 . 

P rn&X. 

"Pn.in 


( 9 ) 


where from Eq. (8) the maximum and minimum values occur at the 
corresponding maximum and minimum values o f cos ^ \ + (j>^ . Thus, 


L z = 


I + R 
I -R 


R = 


L ~ I 
L+l 


( 10 ) 


The phase angle <p follows immediately from Eq. (8) by measuring 
the distance di between the specimen at x=0 and the first pressure 
minimum (see Fig. 2), 


- . - - 1 + 4> = IT 

i-M 0 2 


4> = IT _ 2. Kd t 

l-M/ 


( 11 ) 
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2 . 2 Impedance 


In the absence of a reflected wave (see Fig. 1), the effec- 
tive impedance of a specimen located at x=0 is defined as the ratio 
of the acoustic pressure to axial acoustic velocity. Using this 
definition and substituting Eqs . (4) -(7) yields 


u (o, -t) 



I + Re_^ 
I - Re i4> 


( 12 ) 


Separating Z into its real and imaginary parts, 


g 


l-R z 

+ i 

2 R sin <b 

1 + R 2 - 2 R cos §> 


I + R 2 - 2 R cos $ 




(13) 


resistance 


reactance 


To verify the validity of Eq. (13) , three different check cases 
will be applied. They are (1) rigid wall impedance, R=l,4>=0, (2) 
pressure release surface, R=1 , 4* =TTand (3) no reflected wave case 
(i.e., P 0 c 0 effective impedance), R=4>=0. 


Case 

(1): 

2 

(R-i 

,<J> -o) = '~ r2 - 

1 + R . 

CO 

(14a) 






l-R 


Case 

(2) 

2 

(R = 

l <j> = Tr)= l"R 2 

1 " ^ - n 

(14b) 





(l -t- r) 2 

l -hR 


Case 

(3) : 

2 

(r = 

II 

ll 

•0 

0 


(14c) 


It is clear from the above that Eq. (13) appears reasonable. 

Solving for the resistance and reactance in terms of the 
parameters L 2 and di results in 
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Equations (15) and (16) show that the mean flow reduces the effective 

wavelength of the incident sound via the cos ( — — ' ■■■ term. This 

^ — v I— M 2, / 


can be seen by replacing K = 


Kd, \ _ 


= cos 


Z - rrd 


(l-M a ) 


V— , thus 


= cos 


2. Trd , 

A eff 


where /\e-f-P = ^\ ( I — M 2 ") 


VORTICAL MEAN FLOW CASE - BASIC EQUATIONS 


In a real duct flow environment the duct flow speed is not 
uniform but instead contains the boundary layer region wherein the 
velocity profile varies from zero at the duct walls to an almost 
maximum value at the layer outer edge. The principal effect of 
the boundary layer region is to refract the sound pressure into 
the wall for sound propagating in the same direction as the duct 
flow (the downstream case - see Fig. 3a) and into the duct center 
for sound propagating in the direction opposite to the duct flow 
(the upstream case - see Fig. 3b). 

The basic equations derived in Section 2 governing the propaga- 
tion of sound in rigid-walled rectangular ducts containing a mean 
flow are generalized below to apply to vortical mean flow. Figure 4 
is a schematic of the rectangular duct geometry and coordinate 
system used in the derivation. 
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The equations governing the conservation of mass and momen- 
tum flux in the y and z directions are identical to Eqs. (1), (2b) 
and (2c) respectively. The x-component of the momentum flux con- 
servation equation, however, takes on the following form 


Ail + v G Ail + f A^\ v + w + _L = o 

bt bx \byj \bz J Po bx 


(18) 


The convected wave equation is solved in the usual way to 


become 


^ ^ v J z* df de aTat r °" 


<^M 0 bv bM. dw 


Idy 


bx bz bx 


= O 


(19) 


unrfornA -flow part shear flow part 

The following wave-like solution to Eq. (19) is proposed based 
on a generalization of Pridmore-Brown' s approach, 


r \ 

| / . \ I 

t-XjL) 



P ( x > y. + ><K> V o) 

1 CO V 

= e 


F (y, z-, X-, V 0 ) 

(20) 

v (x,y, a .,t 5 <K;V 0 ) 

ICO 

= e 

(*-§■) 

3(y^>^;Vo) 

(21) 

w (*> y> z i * i K > v„) 

= e iai 

(>■« 

h (y> z >K) v o) 

(22) 


Here < is the system eigenvalue. It will be shown later to 
be related to the group velocity of the propagating sound field. 

Note that the acoustic pressure and transverse velocities are written* 
to depend explicitly upon the mean flow V 0 . The functions g and h 
can be immediately related to F by substitution into the transvers 
momentum equations [see Eqs. (2b)] to yield 

= -*t— 

i (O o co (l - M„) 



) 

— (23) 


10 



( 24 ) 


_ ^ F ^y> > V#') 

h(y,z; V d ) = bz. 

ip 0 co (l-^ M„) 


Substitution of Eqs. (20) thru (24) into Eq. (20) results in the 
following form of the wave equation, 


S. 2. 

dF dF 
+ — : + 


ZK 


by 1 bz (l-^Mo) 


d M 0 dF bM 0 dF 

dy by dz. bz. 




(l-tf M $-#*- 


F=0(9) (25) 


It is convenient to non-dimensionalize Eq. (25) by letting 

F (y;^3 A)V 0 ) = PreF f (S,yiX } V 0 ) } $ =± . J = _E_ (26) 

H H 


where P re £ is arbitrary since Eq. (25) is linear. The quantity 2H 
is the duct height (see Fig. 4). Substituting Eq. (26) into Eq. (25) 
the nondimensional wave equation is written 


d*-P 

bj z + bf 


0-KM.)| 


^kYi£) /jjkVi£ 

\b§Aby + 


y/^dy/J 



(27) 


where KHco/c 0 . 

4. PERTURBATION SOLUTION 

A perturbation scheme to solve Eq. (27) is sought by expanding 
the eigenfunction f and eigenvalue K in a power series of M c , the 
mean flow center-line Mach number as follows 
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f (S , 7 . ■ > V.) = f . (3 , 7 iff) + Me f , £ , 7 , *) t M e * i t Ci , 7 . 4) 


Me f . (j , y> <K) ( 28 ) 


M c 2. ■+■ ' * ' 


M c tfj (29) 


Here the mean flow center line Mach number is defined as 

, 7 = l) = _&_(/•”,*.«) (30) 

The physical meaning of Eqs. (28) and (29) is straight-forward. 
The expansion procedure assumes that the effects of convection and 
refraction are small - sufficiently small so that the leading or 
dominant term is independent of the mean flow. This is the physical 
basis for the perturbation approach adopted. 

4 . 1 Mean Flow Velocity Profile and Boundary Conditions 

Before proceeding further, the mean flow velocity profile must 
be specified. The analysis will be restricted to ducts sufficiently 
long to establish fully-developed mean flow. Thus the mean flow 
boundary layer will be independent of the duct axial position. 

This means that the boundary- layer thicknesses in both the y and z 
direction will extend to the duct center defined in Fig. 4 as 
z=H and y=W. The rectangular duct geometry is such that the mean 



; = o 


<?\ (y°) - ■+• 


Oo 



l - o 
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flow is symmetrical in the y and z direction with respect to the 
duct center. With this understanding, the solution to Eqs. (28) 
and (29) will be derived only for the quadrant defined by z^H and 
y SW. 

The following fully-developed mean flow is proposed. 



(31) 


where N is an arbitrary parameter that depends upon the flow Rey- 
nolds numbers (V c H/£_) and/or (V c W/j^ , and V c is the mean flow 
centerline value. For Reynolds numbers of the order of several 
million, representative of the usual kind of duct turbulence, N=l/7. 
In nondimens ional terms, Eq. ( 31 ) is written 

M (^?) = ".(-a-ViY = M c[ 3^y (32) 


where p =H/W is the duct aspect ratio. Substituting Eqs. (28) -(32) 
into Eq. (27) and collecting coefficients of the powers of M c yields 
to order M c 2 , 


M c 


Y-f 


* 3 Z ^7 




+ 


M C J 


^Y, 








a 3 


7 


+ (y 


O'-CK + Y(f 7) N ^.f. -2 


+ 
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r 



i 3 ™ (syf" 


= o 


(33) 


The first term in Eq. (33) is the classical wave equation govern- 
ing the propagation of sound in a rectangular duct for the ease of 
zero grazing fZow. The solution is well-known and corresponds 
physically to the propagation of all the cut-on modes consistent 
with the source excitation frequency. Observe that the influence 
of the lower -order solutions act as nonhomogeneous source terms 
on the r ight - hand - s ide of the higher-ordered equations. 

The governing sound propagation equations described by Eq. (33) 
must satisfy the rigid-wall boundary conditions of zero normal flow 
defined as, 

° > Z -> A =) = v (x f z,-t) = w(x, y, o,-t ) = w (x,y> H^jaO (34) 


From Eqs.(23) and (24) this translates into the Neumann boundary 
conditions 


£ F (o, z.) 


£F(W,z.) = ZF(y,o) 
dy £ z. 


if(y> H ) = 0 

<)z. 


(35) 
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For the special case wherein the zero -ordered perturbation 
represents the plane-wave sound field, the leading terms of 
Eqs. (28) and (29) become 

($> f) = tfo = 1 (36) 


In this case, Eq. (33) simplifies to 


^ + . 2 ('kh) 1 




(p$?y 


* K. 


(37) 


and 




- — 2 N 


f(£f 


M “ l 




7 J 


+• 


(KH)' 


iT +o ~(p$ f ZN +(V+ -2*,) 


(38) 


4.2 General Solution 


Equations (37) and (38) are the well known Poisson Equations 
subject to the Neumann boundary conditions defined by Eq. (35) 
Solutions to the functions f ( (£> 7 ) and f 2 (J>7) are sought by 
expanding them in a double Fourier series in terms of the homogeneous 
eigenfunctions cos (n?7T^£) cos(£ 7 m^). The general solution to Eq. 

(37) will be formulated. First expand f ( (£,- 7 ) as follows 
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oo oo 


= ^ y c m < h * js $) =« (<”>)) = 


m=0 ^ =0 


Oo 

E 


Am Cos 


m= I 

oo oo 


( m7T / 3 i) + b 2 cos 7* 7) 

■Ul 


( 39 ) 


+ 


^ ^ C mi cos(mvp^ cos (Xny) 


™ = I -dsl 


Next, expand the nonhomogeneous terms (RHS) of Eq. (37) in a similar 
manner 


2. h) (j3 § +■ ^ <K' | 4- 2 (3, 






<q CoS 
1 m 


ms| 


00 

( mTr p$) + 2 (KB) 2 " COS ^ 71 ^ 


i = I 


+ 


2.(KH) 2 


OO 00 



Im^ COS 


1=1 / - 1 


C m ^ z 3 S) 005 (i ^ 


(40) 


The coefficients Ago, Aju, Bj> , Cj^£, Gqq, G m , and Ijh£ are 
constants. The constants G 00 > G m , % and I m £ are defined using 
standard Fourier analysis as 


£ 


H Z 3 / / (/’i?)" d 7 = 

t J (n*>) 


(41) 


' * 


3 m = 2 


Z 3 f J (p$y) Ncos ( ,r 'Y$) e *$ <i 7 = 


- 


ah-i 


0 0 


(n+i) 


(42) 
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where 



0 


( 43 ) 


In a similar manner, the coefficient may be evaluated as 



(44) 


where 



( 



N 

1 


cos 



(45) 


The final coefficient I m ^ is 

= 4 


* ' 


m - 




o O 


The solution to Eq. (37) follows immediately upon connecting the 
coefficients A 00 , A ra , B£ , C m £ to the known coefficients Goo, G m , Pty , 
I m £* Substitution of Eq. (39) into Eq. (37) and collecting coeffi- 
cients of the various cosine terms yields the following relations 


= - 2. ('KH) 2 -g nri , _ -2(Kh') 3 ‘ Hje 

m z TT z p 3 ' 




i Z l T Z 


t ~ 


•7T 2 (m Z A Z + JL 2 -') 


(47) 



The coefficient G 00 is determined by satisfying the compatibility 
relationship defined as 



Thus 



(48b) 


Putting it all together, the first term perturbation solution is 


f (S > ?) = A 


(N+|) 7 T 2 - 



Cos(mTTp§) 
m 2- 


-^Ckh) 2- V ’ 

(N- h f) T r x 

1 = I 




(49) 


The coefficient A 00 in Eq. (49) is arbitrary for the Neumann problem. 
The following tentative physical interpretation is offered. A 
nonzero value couldbe interpreted as a source of acoustic energy. 

The idea here is based on Morfey 1 s le study of production of acoustic 
energy in a flow duct/sound system from interaction between the 
sound field and the mean flow vorticity (JSV: Vol 14, 159-170, 1971). 
Morfey showed that sound can be generated (or absorbed) by the 
mean flow providing the mean flow is (1) vortical and (2) three- 
dimensional. Both conditions are satisfed in the present case. 

Sinoe it is beyond the soope of this study to consider the production 
of acoustic energy , it will be assumed herein that the sole effect 
of the mean shear flow is to redistribute acoustic pressure . Thus 
it follows that the plane wave pressure, normalized to unity, must 
satisfy the constraint that 
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o 


(50) 


/ _/ p(S>7) d i d 7 = 1 J J - f <($’7) d 5 d 7 = 


for all i. 


5. SOLUTION FOR THE CASE N=1 

A detailed solution to second order for the case N=1 is presented 
for a linear mean flow velocity profile. The advantage of the linear 
velocity gradient case (N=l) is that Eqs. (43), (44) and (45) can be 
easily integrated to yield analytic expressions for the perturbation 
solutions of f, (J" , y) and For this flow, the refraction 

effect of the shear flow (e. g . ^ N ) is larger for this case than 

for, say, the N=l/7 velocity profile. To further simplify the analy- 
sis, the duct cross-section will be square so that B=H/W=1 (this 
shape corresponds to the Langley flow duct facility). 

5.1 First-Order Solution 

Substituting N=1 into Eqs. (43), (44) and (45) and evaluating 
the integrals yields 


rri 


(m-rr): 




- 2 . 
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Substituting these values into Eq. (49) results in the following 
solution to f, (£,>?), 



Equation (54) is exact. It is clear that it satisfies the 
boundary conditions 


if = jf, (‘> 7 ) = jf, 6)°) . _ 

^7 ^7 


In order to satisfy the constraint imposed by Eq. (50) that the 
mean flow vorticity does not act as a source (or sink) of acoustic 
energy production, it is straight-forward to show that A oo =0 in 
Eq. (54). With A oo =0, Eq. (54) is shown graphically in Fig. 5.. 

The symbols shown in Figure 5 represent the solution Eq. (54) 
calculated to within one percent accuracy. The dashed lines 
represent the following curve fit, 

(4,7) “ (^h) 2 I . 02.1 + .042 cos cos 2 (¥) 

— (-33 9 + .o^ & ^(. 339 -t-.Oi&^sin 2, s i n z (jU^ 
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As an independent check on the accuracy of the above curve fit 
t0 f i (J;?)* i1: satisfies the energy constraint Eq. ( 50 ) , to the 
following accuracy 

1 1 

J A ($>i) 6 5> ir } = s - 4 * l0 ' s C57) 

O O 

The purpose of the above curve fit is to simplify enormously the 
calculations required to compute the second-ordered perturbation 
solution The curve fit is sufficiently accurate to re- 

present both the function f, (£,>7) and its first derivat ives 
and '/^7 both of which are required to calculate f a (£,^) as shown 
in Eq. (38). It is clear from Fig. 5 that the effect of the shear 
flow is to redistribute the initially plane wave acoustic pressure 
in such a way as to increase (decrease) the local pressure near the 
wall for downstream (upstream) sound propagation (i.e., near ^ =0 
and rj = o ) relative to its value at the duct center (J = ^_/). 
From Eq. (48) the eigenvalue 





(58) 


Further physical interpretation 


(s > y) i 


is obtained, 


is deferred until the solution to 


5 .2 Second-Order Solution 


The solution to the second-order pressure f z (^ ; ^) perturbation 
is similar to that used to calculate f, (5,77). Substituting N=/3=l 
into Eq. (38) results in the following Poisson equation 




+ 




= -*5 + ( kh ) : 


if, 

drj 






( 59 ) 
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The eigenvalue <K Z is completed from the RHS of Eq. (59) 
posing the compataiblity relation that 

I 

(RHS) dy = O 

o 0 



by im- 


(60) 


Substituting Eq. (56) into Eq. (59) and integrating results in the 
following approximate expression for , 


i\ z ^ -.015 -v- .009 


(61) 


Thus an approximate solution to the effects of refraction on the 
axial propagation of the plane wave, valid to second-order, is 


<K - 



M c 


0\5 — . oo 9 




(62) 


Referring back to Eqs. (4) -(6), the eigenvalue i\ is related to 
the propagation speed of the sound energy - that is, the group speed 
(Vg) . Recall that for sound propagation through a dispersive media, 
the group speed is a measure of the average rate of transport of 
acoustic energy. Its mathematical definition is 


where 


V« = 


_ Z 


60 


3K 
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I - M c z 
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.0I5-.009(KH) Z 
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Performing the necessary differentiation, the group speed, normalized 
to the local sound speed c 0 , is approximately 



(65) 


Figure 6 shows a comparison between the group velocity defined 
by Eq. (65) and the group velocity corresponding to the propagation 
of -plane wave sound in a uniform mean flow. It is well known that for 
sound propagation in a uniform mean flow, the group velocity is 
equivalent to the sound wave phase speed* Thus 

-4- = = 1 +M c (66) 

*-•0 C 0 


Figure 6 shows that the group speed for the- shear flow case is 
quite insensitive to the centerline mean flow speed 
in contrast to the uniform flow case. The effects of refraction 
appear to focus the group velocity to values near c Q for both up- 
stream and downstream sound propagation. 

With <K Z specified, the RHS of Eq. (59) is known. Its solution 
follows exactly the procedure described above in Section 5.1. 
Substituting Eq. (56) for f ( ($, 7 ) into the RHS of Eq. (59) results 
in the following rather complicated Poisson Equation, 
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= (^Kh) z J 7r y sin (^r £) 
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4 ^.339 -+• .016 y) Yj 
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(.33 9 + . 016 ,§)('339 4 • < 3(6 >7) sin 2 


.033 -_L 
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(£7)-(s>?y 


+ (k.h) + J .018 t( 2 $^-±) 


. 02.1 4 .042 cos 2 


c “^) 



- ^.339 4 . Oi ^ 339 4 . 016 rT) sm Z f_ n $\ sin 2 / 7r 

' \ 2. / V ^ 

+ (KH)*. B(S, 7 ) 




(67) 


Obviously, the RHS of Eq. (67) is quite complicated. To reduce 
the computational calculations, it has been simplified as shown in 
Figures 7 and 8 to the following curve fit expressions for A(£,^) 
and B(£,^) , 


A (£> j ) - - 033 + ( Sy ) -^73 - 1.533 (£7) -.5-67 


( 68 ) 


B (§>j) - --OI 5 .09 7 -.236 ( 5 ^ + -021 (£ + y) 


(69) 


As a check on the accuracy of Eqs. (68) and (69) , the compatibility 

relationship requires that the integral of the RHS of Eq. (67) 

vanish yielding 
1 1 


(K H ) l A(i-7) + ( l< H)' f B(i. 7 )|dSd 7 = .oi(KH)\003(kH) 4 (70) 


o o 


24 



The RHS of Eq. (70) is sufficiently close to zero to verify the 
accuracy of the curve fits defined by Eqs. (68) and (69) . Using 
these equations and following the approach used in Section 5.1 above 
to evaluate the Fourier coefficients to the first perturbation pres- 
sure solution, the solution to the second perturbation pressure 
£ 2 (5>7) is 
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Equation (71) is rewritten in the more convenient form with one 
term proportional to (KH) 2 and a second proportional to (KH) 4 , 


=(kh) 2 c($,7) + (kh)* 


(72) 


The coefficients and approximate curve fits to c(g,>y) and 
D(£’7^ are shown graphically in Figs. 9 and 10. The symbols represent 
numerical solutions to Eq. (71) and the dashed and continuous lines 
represent the following curve fits, 

C fe>7)“ ^(71^) cos(TJrj') + .OlQ ( 73a ) 


^ 7 ) - ,001 c °s^7r^ cos(V^ -t-.ooiS sin sin (ft (73b) 

- - oz (£ -00^3 (\-£)(\-rj) ~ .oo\Q (\-^y^\-rjY 

Equations (73a) and (73b) satisfy the constraint of zero energy 
transfer between the mean vortical flow and the sound field. The 
integrals of the curve fits, Eqs. (72) and (73) , are 


* I 

J [ - f z($-’j) d $ A 7 

0 o 


.001 (kh) Z - ■ 0008 (kh) 4 


(74) 


which approximate zero sufficiently accurately. Thus the free 
constant equivalent to A 00 in Section .1 above is set equal to zero. 

Combining Eqs. (28) , (36) , (56) , (72) , (73) and (74) results 
in the following approximate solution to second-order, for the sound 
pressure 
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(75) 


Here p ref is an arbitrary source strength. Equation (75) is only 
approximately correct - but is is sufficiently accurate to 
demonstrate the importance of the effects of refraction by the mean 
flow. Figure 11 shows the acoustic pressure distribution across the 
duct for the case KH=1 and Mc=0.2. The refraction effects are quite 
small. Figure 11 suggests that the ratio of the acoustic pressure 
at the origin y=z = 0 to its value at the duct center y= z=H(£ = y = ! ) 
is a good measure of the distortion of the plane wave due to the mean 
flow refraction. From Eq. (75) , 


1 ■+■ M c 

.0 63 (KM) 2- 

+ Me 1 

.02+(kH) z + .003S ("kh) + 


1 + M C 

-.I05(kH) j ' 

] + M c 
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" .018 (KH) 2 + . 0007 
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Figure 12 summarizes Eq. (76) for a wide range of mean flow 
center Mach numbers and sound frequencies. It is clear that for 
low speed flow and for sound frequencies below the first mode cut- 
. on frequency (corresponding to a axial wave number value KH< w /2), the 
effects of mean flow refraction are on the order of ldB or less. 

For application to the Langley SWT duct flow facility, clearly 
the effects of refraction are unimportant. For downstream sound 
propagation, Fig. 17 shows the effects of the mean flow is to re- 
fract the sound into the duct walls. For the upstream case, the 
sound is refracted into the duct center. 

It is implicit in the formulation of the perturbation scheme 
described by Eqs. (28) and (29) , that both M c and KH must be small 
in some "sense". The frequency parameter KH is straight-forward. 

For perturbations about plane-waves, KH<?r/2. The range of validity 
of the centerline Mach number is defined in terms of convergence 
of Eqs. (28) and (29). From Eq. (78), it is clear that the maximum 
value of Me for which Eq. (78) is valid is related directly to the 
value of (KH) . Thus, Me can be larger for smaller values of (KH) . 
Consider for example, the numerator of Eq. (78) . For P QO to be 
valid, the following two constraints must be met, 

• 0&3 (kh) M c K \ (79) 



. 024 -t- .0035 (kh) 2 


< .063 (kh) 2 Me 


(80) 


In general, the m^+i term must be smaller than the m^^ term. Equations 
(79) and (80) show that for KH<7 t/ 2, all subsonic center duct flow 
Mach number Values are permissible. 
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6 . 


CONCLUDING REMARKS 


The perturbation model shows clearly that for low speed 
mean flows and for sound frequencies corresponding to plane -wave 
excitation, the mean flow velocity gradients only negligible 
distort the plane wave transverse pressure distribution across 
the duct. This conclusion is conservative in the sense that 
for mean velocity profiles defined by Eq. (32) , the velocity 
gradient (vorticity) in either the y or z direction is proportion- 
al to N; thus selecting N=1 results in a much larger value of 
vorticity than the usual turbulent value of N=l/7. It is im- 
portant to remind the reader that the above conclusion has been 
verified only for perturbation about the -plane-wave ease. An 
important finding of this study is that the effects of refraction 
do not significantly affect the standing wave pattern in a 
wave guide. This is valid providing the flow speed is low and 
only plane wave motion is excited within the waveguide. 

A rather interesting result of this study is that the 
mean flow velocity gradients seems to focus the propagation 
of acoustic energy into a narrow window around the undisturbed 
adiabatic sound speed, both for upstream and downstream propaga- 
tion. This is in contrast to the uniform mean flow sound energy 
propagation speed of c 0 (l+M 0 ), which is a strong function of the 
mean flow velocity. 

Another interesting result is that the perturbation tech- 
niques permits the possibility of allowing for the production 
(or absorption) of acoustic energy through interaction between 
the mean (vortical) flow and the sound field. To consider this, 
however, the energy conservation equation must be studied using 
the technique described by Morfey. 
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Figure 1. Schematic o£ SWT Impedance Model 
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FIGURE 3b. 


SCHEMATIC OF UPSTREAM SOUND PROPAGATION 





Symbols : Numerical Calculations 'v V N ^ 

A^roxima-te Curve. Ftt: 

-.10-- ft (i> 7) - - 042 cos 2, co,s2 ( 7r -^-) 

-(.339 +.016^.339-1-. 016 ^)sin a (^4.)x 

Sin a ( 7 f ^^-+- .021 

FIGURE 5. CURVE FIT OF FIRST PERTURBED ACOUSTIC PRESSURE 
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Mean Flow Velocity Profile. 
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FIGURE 7. CURVE-FIT TO THE TERM PROPORTIONAL TO (KH) 2 ON THE RHS 
OF EQ. (67) 











FIGURE 9. 


CURVE-FIT TO THE TERM PROPORTIONAL TO (KH) 2 ON THE RHS 
OF EQ. (71) 
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ooi <x3s(7r£)cos(7ry)+ .0015 sin(n^) sin (V yj) - . 02 £ rj (j.2 

+ - 0043 (|- j )(|- y) T . 0 018 (' I -^) 2 (|-^) 2 ' 


FIGURE 10. CURVE-FIT TO THE TERM PROPORTIONAL TO (KH) 4 ON THE RHS 
OF EQ. (71) ... v.... 
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FIGURE 11. EFFECT OF REFRACTION ON THE ACOUSTIC PRESSURE DISTRIBUTION 
ACROSS DUCT FOR THE CASE KH=1 AND M =0.2 






FIGURE 12. EFFECT OF FLOW AND SOUND FREQUENCY ON THE RATIO OF 
WALL TO CENTERLINE ACOUSTIC PRESSURE 
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